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I Abstract. We present 3D hydrodynamical simulations of the separated bi- 

I , nary RS Ophiuchi (RS Oph), a recurrent nova and potential progenitor of a SN 

O ' la. RS Oph is composed of a red giant (RG) and a white dwarf (WD) whose mass 

|H I is close to the Chandrasekhar limit. In an isothermal scenario, the WD accrets 

c/3 . about 10% of a 20 km/s RG wind by a non-Keplerian accretion disk with strong 

spiral shocks, and about 2% of a 60 km/s RG wind by what we term a 'turbulent 
accretion ball'. A significantly larger impact have the thermodynamics. In an 
adiabatic scenario only about 0.7% of the 20 km/s RG wind is accreted. The 
^ ' rate of change of the system separation due to mass and angular momentum loss 

OO \ out of the system is negative in all three cases studied, but is ten times smaller 

0^ ■ for a fast RG wind (60 km/s) than for a slow RG wind (20 km/s). The results 

\l ' demonstrate that existing nova models and observed recurrence times fit well 

together with 3D wind accretion and that RS Oph is one of the most promising 
i systems to become an SN la. 

o 



^ ' 1. Introduction 



Type la supernovae (SNe la) are cornerstones of modern cosmology as a measure 
for cosmological distances, and they are crucial building blocks of the universe, 
as production sites of a large part of iron group elements. In an SN la a white 
dwarf (WD) explodes after his mass has surpassed the Chandrasekhar stability 
limit. The existence of SN la thus is closely linked to the question of mass gain 
of the WD by accretion. The layer of accreted matter occasionally undergoes 
a thermonuclear runaway (nova) and much o f the previously acc reted matter is 
blown into space again. Based on ID models, lYaron et all (|2005l ) find that a net 



mass gain of an already massive WD {M ~ IAMq) occurs only for accretion 
rates above roughly lO~*M0yr~^. Moreover, to achieve a net gain of 0.05 Mq 
the WD must undergo many nova cycles, requiring a total time on the order of 
1 to 10 million years. 
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P articularly intere sting binary star systems in this context are recurrent 
novae ( Anupamal |2002| ) . From the short decay time of their nova light curves 
the existence of an already massive WD close to IAMq can be inferred. In the 
subclass of symbiotic-like recurrent novae, which we are interested in here, the 
secondary is a red giant (RG) star which undergoes a substantial mass loss in the 
form of a stellar wind. The RG typically does not fill its Roche lobe and mass 
transfer onto the WD occurs through wind accretion. The accretion efficiency 
then is on the order of 10% of the RG mass loss. Nevertheless, with a physically 
plausible RG mass loss rate of 10~^MQyr~^ the necessary WD accretion rates 
can be reached. Sustaining these conditions over several million years is not in 
a priori contradiction with typical RG life times. Symbiotic-like recurrent nova 
thus are plausible candidates for SN la progenitors, but a number of questions 
remain open. 

The best observed system of this kind is RS Ophiuchi (RS Oph), which we 
investigated already in IWalder et all ( 20081 ) where we presented a 3D hydrody- 
namical simulation of the recurrent nova RS Oph during the accretion phase and 
nova outburst. Two issues we take up here again, in the light of new simulation 
data: accretion rates in a recurrent nova binary star system and the competition 
between mass loss and angular momentum loss out of the system. The latter is 
crucial as it affects the separation of the system and thus the accretion rate. 



2. Model Problem and Numerical Method 



The symbiotic-like binary star system RS Oph consists of a red gia nt and a white 
dw arf. The orbital period is 455 days (Fekel et al. 2000; Dobrzvck a aiid Kenvon 
1994). The system undergoes a nova outburst about every 22 years (jAnupama 
2OO2I ). For the most recent out bursts in 19 85 and 2006, exquisite panchro matic 



observational data ar e available (IShqre et al. 1996; Bode et a/. 2006 : Sokoloski et al 
20061: iDas et a/.ll2006l: [O'Brien et a/.ll2006l : lEvans et al.ll2006l : IWorters et a/.ll2007l : 



Hachisu et a/.ll2007l : TBode et a/.ll2007l ). We adopt masses of 2.3 M© and 1.4 Mq 
for the RG and WD, respectively and a corresponding separation between the 
components of a = 2.68-10^^ cm. We assume phase locking of the RG, a RG ter- 
minal wind velocity of either vug = 20 km/s (simulation slow) or v^q = 60 km/s 
(simulation fast) in the rest frame of the RG, a mass loss rate of 10"'' MQyr"^, 
and a wind temperature of 8000 K. The radius of the RG is smaller than its 
Roche lobe, accretion by the WD occurs from the RG wind. 

We solve the Euler equations in 3D with a nearly isothermal polytropic 
equation of state with 7 = 1.1, resulting in a thermal structure that comes 
close to that obtained from photoionization models of related symbi otic binary 
systems ( Nussbaumer and VogelllT987l : iNussbaumer and Waldeijll993l ). For com- 
parison, we also did one adiabatic simulation with 7 = 5/3 and -^rg = 20 km/s. 

All simulati ons were performed with the hydrodynamical codes from the A- 
MAZE package ( Walder and Folini 2000l ). The accreting system was relaxed to 
a quasi-stationary state. Mass and angular momentum losses out of the system 
were then computed through spherical shells arou nd the center of ma ss. For 
further details on the simulation set up we refer to lWalder eta!\ (l2008l ). 
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Figure 1. Density in orbital plane for slow (left) and fast (right) simulation. 



3. Large scale patterns, accretion, system evolution 

During the roughly 22 years of quiescence between successive nova outbursts, 
the circumstellar environment gets filled with RG material, which ultimately 
reaches out to several 10^^ cm. The matter distribution on these large scales 
differs markedly between the slow and fast simulation, as shown in Fig. [TJ In 
the slow case, the situation is fairly smooth and spherically symmetric. The 
spiral pattern in the equatorial plane is rather faint. A density contrast of a 
factor 2-3 between polar (less dense) and equatorial directions exists. This is 
sufficient to transform a symmetric nova outburst of the WD into an asy mmet- 



ric nova remnant with an aspect ratio of about 2:1 (IWalder et al.l 12008). Th e 



observed nova remnant in RS Oph is asymmetric as well (jO'Brien et al. 



SqqI), 



a quantitative comparison remains to be done. For somewhat simpler underly- 
ing models, observable quan tities have already been modeled and s uccessfully 
compared with observations ( Vavtet et a/.l l2007l : lOrlando et ad 1200 9'). The fast 
simulation shows a much more pronounced, farther reaching spiral structure in 
the equatorial plane (Archimedian spiral). Density contrasts across spiral arms 
exceed one order of magnitude even at distances ten times the system separation. 
Densities in the polar direction are much more smooth. 

The fast and slow simulation also show substantial differences in the ac- 
cretion geometry around the WD, as shown in Fig. [2l In the case of the slow 
RG wind, accretion permanently takes place through a non-Keplerian disk in 
which angular momentum is transported via spiral shocks. In the case of the 
fast RG wind, the vicinity of the WD can rather be characterized as a super- 
sonically turbulent accretion ball. Matter flows in from various directions and is 
redirected, slowed down, and thermalized by a complicated and rapidly chang- 
ing network of shocks, while angular momentum is transported outwards. For 
a more thorough discussion of supersonic turbulence in a non-periodic setting 
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see iFolini and Waldeil ( 20061 ). The overall appearance of the vicinity of the WD 
is much more fluctuating than in the slow case. Disk-like accretion states also 
occur, but only as transient features. In the slow wind case, densities around 
the WD are up to two orders of x nagnitude larger in the orbital plane than 
perpendicular to the orbital plane ( Walder et al. 20081 ). In the fast wind case 
no such pronounced density contrast exists. The accretion rates change from 
around 10% Mrg in the slow case to roughly 2% in the fast case. By contrast, 
the accretion rate is a mere 0.7% in the adiabatic simulation we performed. 

In all the cases considered here, the cumulative effect of the system's mass 
and angular momentum losses points to a shrinking of the binary orbit. The loss 
of angular momentum more than outweighs the mass loss of the system and the 
two components spiral inwards. This picture is dominated by the quiescent phase 
of RS Oph, as only much less matter (around 10%) is involved in the occasional 
nova outburst. Absolute numbers turn out, however, to depend strongly on the 
RG wind velocity. Expressing the cumulated losses in terms of da /a with a the 
system separation, we find da/a ~ 8 • 10~® yr~^ in the slow case but about ten 
times less in the fast case, da/a ^ 8 - 10~^ yr~^. Of much less relevance here are 
the thermodynamics. For v^q = 20 km/s and 7 = 5/3 we find da/a ~ 6 • 10~^. 



4. Discussion and Conclusions 

Our results further support the idea that the recurrent nova RS Oph is a good 
candidate for a progenitor of an SN la. In particular, the mass transfer rate 
required by nova models to reach a net mass gain over several nova cycles can 
be reached for reasonable system parameters. Also, the orbit of RS Oph is found 
to shrink under current conditions. Chances to maintain or even enhance the 
current mass transfer rate in the future thus are intact. 

However, the three critical time-scales which determine the ultimate fate of 
the system are all of the same order. RG evolution models predict a high mass 
loss for a RG for some million years only. This time scale is similar to what is 
needed to shrink the orbit significantly and to the time necessary for the WD 
to accrete enough mass to become unstable. 

To further investigate the possibility of RS Oph finally becoming a SN la 
simulation results should cover the long term evolution of the system and its 
components, over several nova cycles up to several million years. In addition, 
multi-dimensional nova models with better predictions of the mass gain condition 
of the WD would be most welcome. On the observational side, better estimates 
on the ejected mass in a nova outburst would be highly desirable to test the 
consistency with existing nova models and related accretion physics. 

Acknowledgments. The authors wish to thank the crew running the HP 
Superdome at ETH Zurich and the people of the Swiss Center for Scientific 
Computing, CSCS Manno, where the simulations were performed. 
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Figure 2. Accretion around the WD in the orbital plane (left column) and 
perpendicular to the orbital plane (right column) for vrq = 20 (top row) and 
^^RG = 60 (middle and bottom row). Shown is density, logarithmic scale. In 
the vyig = 20 case, occurs permanently through a non-Keplerian disk with 
spiral shocks. In the wrg = 60 case, accretion occurs predominantly through 
a 'turbulent accretion ball' (middle row) but occasionally also through a non- 
Keplerian disk (bottom row), similar to the wrg = 20 case. The spatial scale 
of each of the above panels is roughly 10^^ cm. 



